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Magneto-optical imaging technique for hostile environment: the ghost imaging 

approach. 
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We develop a new approach in magneto-optical imaging (MOI), applying for the first time a 
ghost imaging (GI) protocol to perform Faraday microscopy. MOI is of the utmost importance 
for the investigation of magnetic properties of material samples, through Weiss domains shape, 
dimension and dynamics analysis. Nevertheless, in some extreme conditions such as e. g. cryogenic 
temperatures or high magnetic fields application, there exists a lack of domains images due to the 
difficulty in creating an efficient imaging system in such environments. Here we present an innovative 
MOI technique that separates the imaging optical path from the one illuminating the object. The 
technique is based on thermal light GI and exploits correlations between light beams to retrieve 
the image of magnetic domains. As a proof of principle, the proposed technique is applied to the 
Faraday magneto-optical observation of the remanence domain structure of an yttrium iron garnet 
sample. 


Introduction. The magnetic behavior of a material is 
directly correlated to its domain structure (DS). The DS 
is determined by, amongst other things, chemical compo¬ 
sition, microstructure, growth processes and post-process 
treatments. The capability to observe the DS is funda¬ 
mental for designing tuned material properties and to ob¬ 
tain the desired magnetic behavior [l|. The understand¬ 
ing of the structure and dynamics of magnetic domains 
and domain walls is of the utmost importance in a num¬ 
ber of applications, such as thin-film recording heads in 
the magnetic recording industry and spin-electronic 
devices in information technology [3|. Data about mag¬ 
netic anisotropy, internal stresses, magnetic losses and 
domain wall tWnamics can be inferred by magnetic struc¬ 
ture imaging [3,@. Several imaging techniques have been 
developed over the years, e.g. magneto-optical Kerr and 
Faraday microscopy, magnetic force microscopy (MFM), 
magnetization sensitive electron microscopy, scanning 
Hall probe microscopy, scanning SQUID microscow, and 
Bitter decoration. Magneto-optical techniques [a Q j ex¬ 
ploiting the Kerr effects in reflection and the Faraday 
effect in transmission, in particular, have some advan¬ 
tages, among which the ease of setup and a very good 
time resolution, allowing the investigation of fast dynam¬ 
ical processes for the characterization of magnets and 
superconductors Due to the Faraday effect, the po¬ 
larization plane of the light transmitted through a mag¬ 
netized material endowed with spontaneous local magne¬ 
tization is rotated of an angle @ ~ where L 

is the thickness of the magnetic material, oj is the light 
frequency, c is the speed of light, Mj. is the local magne¬ 
tization component along the wavevector k and e is the 
dielectric permittivity of the medium. The sign of the 
angle (j) depends on the sign of the magnetization. Simi¬ 
larly, due to the magneto-optical Kerr effects, a rotation 
occurs when the li^t is reflected from the surface of a 
magnetized sample [O]. In a Faraday (Kerr) microscope, a 
beam of polarized light is exploited to perform the imag¬ 
ing of adjacent domains, with different magetization ori¬ 
entation, sending the transmitted (reflected) light to a 


spatially resolving detector (e. g. a CCD) preceded by a 
polarization analyzer. Such setup guarantees a magneto¬ 
optical contrast between adjacent domains. The DS, ob¬ 
served by means of the magneto-optical effect, can be 
considered as an object with transmission T(Mfe(x)) (re¬ 
flection R) depending on the local magnetization at the 
coordinate x. In the last years, with the advent of mag¬ 
netocalorics earlier and spincaloritronics then (l^ . the 
possibility to perform magneto optical imaging in tem¬ 
perature variable experiments is acquiring in perspective 
more and more interest. This means, unavoidably, to 
look for the possibility to integrate a magnetic imaging 
system in experimental setups where the optical access to 
the sample is strongly limited, due to temperature con¬ 
trol additional stages [ll| or when the sample would be 
cooled at extremely low temperatures. Until now, MOI 
at temperature down to IK has been achieved by means 
of cryostats endowed with windows, which allow polar¬ 
ized light to reach the sample. However, radiation losses 
due to the presence of windows hinder the efficiency of 
the thermal isolation [H, limiting the performance of 
the cryostat in sub-Kelvin regime. Furthermore, errors 
can be introduced due to the Faraday rotation in the mi¬ 
croscope objective, with visible effects already at applied 
fields of the order of hundreds of mT, which unavoidably 
worsen when approaching the regime of several Teslas 
[l^ . Similarly, the applied magnetic field intensity could 
affect the correct operation of the camera electronics. In 
this paper, we propose an innovative, flexible way of per¬ 
forming MOI of DS exploiting the ghost imaging (GI) 
technique [l3 - [^ . Compared to conventional techniques, 
in a GI setup the imaging system is placed in a spatially 
separated optical path with respect to the one contain¬ 
ing the sample under test; the light-beam that interacts 
with the sample is collected by means of a bucket detec¬ 
tor (a detector without spatial resolution) like, e.g., an 
optical fiber connected to a photodiode. A second beam, 
locally correlated with the probe beam, does not interact 
with the sample and is sent to a spatial resolving detec¬ 
tor (CCD camera). The technique retrieve the image of 
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the DS exploiting correlations [2ll - l^ between the out¬ 
put of the bucket detector with each pixel of the CCD. 
This setup paves the way to realize MOI in extremely 
high magnetic fields, in window-less cryostats or when¬ 
ever the accessible volume in the proximity of the sample 
is limited. 

Theory of thermal Ghost Imaging. Thermal GI (Fig¬ 
ure [1]) requires a couple of spatially incoherent, locally 
correlated thermal beams [Ij, dubbed here as a- 
beam and 6-beam, generally obtained by splitting a sin¬ 
gle pseudo-thermal beam through a beam splitter. In the 
two respective transverse planes, ideally divided in dis¬ 
crete resolution cells (each one larger than a ” spatial co¬ 
herence area”) of coordinate {i = a, 6), the two beams 
present thermal intensity fluctuation statistically inde¬ 
pendent cell-by-cell, 


{{5n^ff) = {n^fj ^1 -I- 




M 


M 


( 1 ) 


where (nx^) and M represent the mean number of pho¬ 
tons and the number of independent spatio-temporal 
modes collected within a cell during the acquisition time, 
respectively. = n^i — (rixi) is the photon number 

fluctuation. The last approximation in Eq. o, valid 
when (uxi) /M ^ 1, is fulfilled in typical experimental 
configurations. The mean photon number is conveniently 
expressed as 


(«xi) = Tj(xi)A i = a,b (2) 


where A is the mean number of photons impinging on the 
resolution cell while Ti{xi) represents the overall locally 
defined transmission/detection efficiency of the channels. 
In GI o-beam first interacts with the object with a struc¬ 
tured transmission profile ~ Ta(xa), then is detected as a 
whole by the bucket detector, i.e. a detector without spa¬ 
tial resolution, whose output signal is Na = 
sum extends over a defined number of resolution cells). 
The 6-beam is detected by a spatially resolving detector 
at the correlated transverse plane, like a 2D-pixels array. 
The output signal from each pixel in the position x;, is 
denoted by The point-by-point spatial correlation, 
that emerges from the splitting of the intensity fluctua¬ 
tion, can be expressed in terms of the covariance of the 
detected photon number as: 


= T6r,(x',)^6x'„x, 

( 3 ) 

where we have considered a nniform transmission- 
detection efficiency in the channel 6, i.e. Tb(x6) = Tf,. 
Therefore, the covariance of a single pixel in Xf, with the 
signal from the bucket detector in channel a turns out to 
be 

(diVadrix,,) = (^Ux'^du-xt) = TbTa(x(,) —, (4) 


which permits to reconstruct on channel 6 the transmis¬ 
sion profile of the object interacting with beam a. This 



FIG. 1. Thermal Ghost imaging schematic representation: 
two correlated beams, a-beam and 6-beam, are generated by 
sending an incoherent beam through a 50% beam splitter 
(BS). Then, beams are sent to two distinct optical path, one 
containing the sample to be imaged and a bucket detector, 
the other one containing a spatial resolving detector. The 
image of the sample is retrieved correlating the output of the 
two detectors. 


is the essence of the GI technique. Eq.s m and O show 
how GI reconstruction emerges point-by-point as a corre¬ 
lation of pairs of resolution cells immersed in a large num¬ 
ber of uncorrelated contributions. Thus, the statistical 
ensemble (namely the number of acquired photons rix^ 
and of values Na) should be ’Targe enough” to allow the 
emergence of such correlation from the uncorrelated con¬ 
tribution; this increase the total duration of the experi¬ 
ment. Obviously, it is of utmost practical importance to 
quantify what does ” large enough” mean for a specific GI 
application. As usual, the theoretical mean values are es¬ 
timated by an arithmetic average of K samples from the 
same population, i.e. {X) —>• E[X] = A^. In 

particular, the covariance in Eq. (jd]) is estimated by the 
statistical quantity Gx^ = E[{Na-E[Na]){n:^^-E[n:^^])]. 
As expected, its value converges to the theoretical one, 

(Gx,) {6Nadn^,) (5) 

with fluctuations around the mean value of {K >> 1) 

(((5GxJ') ~ K-\{6Na6n^,)^) - (SNaSn^,)^ (6) 

«A-i((<5iV„)2)((fexJ^) 


The approximation in the second line of Eq. ([6]) holds 
since almost all the light collected by the bucket detector 
is uncorrelated with the pixel of the spatially resolving 
detector. As expected, the fluctuation scales towards zero 
when increasing the sample size K. Let us consider now 
an object defined by two regions Saj {j = +, —) charac¬ 
terized by two different transmittance Taj for Xq S Saj- 
On channel 6 we can identify two regions Sbj, each one lo¬ 
cally correlated with the corresponding Saj- The signal- 
to-noise ratio (SNR) (sometimes referred as to ’’contrast- 
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FIG. 2. (a) The magneto-optical GI experimental setup for 
Faraday microscopy of magnetic domains. In the scheme BS 
= 50-50 beam splitter, PR. PI and P2 are polarizers, L, LM 
and LI are lenses of focal length / = 100 mm, /m = 100 
mm and fi = 75 mm respectively, (b) The auto-correlation 
measurement for evaluating the coherence area, before and 
after the microscope with magnification factor M = 0.2. 


to-noise ratio” in GI literature 


defined as 


SNR = 


\{C+-C.) 


y'{isc+r) + {{sc.y 


(7) 


should be larger than 1. In fact, by using Eq. ([5]) with the 
substitution of Eq. (HJ one gets {Cj) = /M. In 

the same way, the result in Eq. (j6|) with the substitution 
of Eq.s ([T]) and (dJ returns the evaluation of each noise 
component at the denominator of the SNR as {{SCj^) ~ 
K- 


leading to: 
SNR Ri Vk , 


Tl+ 


R+ + T^_R- 


( 8 ) 


SNR only depends on the transmittance of the two re¬ 
gions Taj, on the number of resolution cells Rj re¬ 
spectively transmitted through the object and on the 
size of the statistical sample K. In particular, in the 
magneto-optical application presented in this work, we 
have roughly = i?_ = R. In this case it is easy 
to show that the SNR is maximized where one of the 
transmittance is null, reaching the optimal value SNR = 

VWi- 

Faraday ghost imaging microscopy: experimental 
setup. The experimental setup of our magneto-optical 
GI experiment is depicted in Figure [2] (a). The two lo¬ 
cally correlated thermal beams, a-beam and 6-beam, are 
obtained by means of a single pseudo-thermal beam di¬ 
vided by a 50-50 beam splitter. Polarized pseudo-thermal 
statistics light is engineered by sending coherent light 
pulses, such as a polarized pulsed laser beam, on a scat¬ 
tering medium with random scattering centers distribu¬ 
tion, in order to obtain a speckle pattern with (pseudo) 


thermal intensity fluctuations. In our case the coher¬ 
ent source is the second harmonic of a Q-switched Nd- 
YAG laser (532 nm), with pulse duration of 10 ns and 
repetition rate of 12.4 Hz. A spatial filter removes non- 
Gaussian components of the laser. Pulses are sent to a 
rotating ground-glass disk, known as Arecchi’s disk [2^ 
(with grain dimensions of some fj,ms). The rotation speed 
is set in order to be considered frozen within the single 
pulse duration, but so that different pulses experience dif¬ 
ferent independent grain distribution. Thus, each pulse 
generates an independent intensity spatial distribution, 
uncorrelated with the previous one. A far-held speckle 
pattern is obtained with lens L, with / = 100 mm, placed 
in / — / conhguration with respect to the disk. A non¬ 
polarizing beam splitter separates the pseudo-thermal 
beam in two balanced, polarized and correlated pseudo- 
thermal beams, with spatial correlation of 98% [^. The 
a-beam is the transmitted beam while the 6-beam is the 
rehected one. a-beam is further polarized by means of 
the polarizer PI and is sent to the sample, a Him of YIG 
of chemical composition (YSmCa) 3 (FeGe )5 0 i 2 , charac¬ 
terized by out-of-plane anisotropy and an average domain 
size of 160 ^m [^. The polarizer P2 allows magneto¬ 
optical contrast observation between adjacent domains. 
In our sample, which is a simple case of uniaxial sys¬ 
tem at the remanence, the magnetization points along 
the direction perpendicular to the sample plane. Adja¬ 
cent domains have opposite magnetization, hence trans¬ 
mittance T, after proper polarization selection, can be 
considered as T{(j)+) = T+ or T{(j)-) = T_ in the case of 
positive or negative Faraday rotation angle respectively. 
In our setup, both a and 6 beams are sent, by means of 
imaging lenses LI of // = 75 mm, to different portion of 
a GGD, with frame acquisition synchronized with laser 
pulses. Bucket detector condition is fulHlled by summing 
GGD pixels in the region illuminated by a-beam. The 
camera is an EMCCD (Andor Luca R camera), with high 
sensitivity and active area of 8 x 8 mm^ and squared pixels 
of size 8/rm. For each laser shot, the output of the sys¬ 
tem is the number of photons Na detected by the bucket 
and the matrix of detected photons corresponding 
to the pixel counts in the region of the GGD illuminated 
by 6-beam. The speckle dimension corresponds to the 
resolution cell dimension, and hence, to the resolution of 
the reconstructed image. The diameter of the speckles is 
evaluated as the full width at half maximum (FWHM) 
of the auto-correlation coefficient estimated as : 


=(«) = Y 


\/[(5nxJ2[,5nx6-H{]2 


(9) 


where ^ is the two-dimensional shift of the selected re¬ 
gion with respect to itself and Xh is the vector position 
of the pixel of the region 6 (see Figure [2] (b)). In order 
to enhance the resolution, the speckle dimension was re¬ 
duced by increasing the diameter of the collimated coher¬ 
ent source impinging on the Arecchi’s disk (see Fig- 
ure[2](b)). Also, a microscope with lenses LM (/m = 100 
nm) able to further reduce the speckle diameter with a 
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a) b) 



FIG. 3. (a) Ghost imaging of magnetic domains of a portion 
(a (400 X 560) fim^ area) of our YIG sample in comparison to 
(b) the direct Faraday imaging of the same portion 


magnification factor M = 0.2 was constructed. In this 
way the number of speckles on the sample was increased. 
The setup finally provides a spatial resolution of 24 ^m 
at the plane of the object, largely sufficient to observe 
the magnetic domains of our sample. 

Ghost imaging of magnetic domains. The ghost image 
is reconstructed evaluating the experimental covariance 
C(x;,) averaging on the total number of frames K ac¬ 
quired by the CCD. We set polarizer P2 angle, minimiz¬ 
ing transmittance r_. In this way we are able to get as 
close as possible to the ideal condition for the SNR, as 
discussed after Eq. We finally set Tf, with a polarizer 
PR in order to avoid saturation of the CCD sensor; any¬ 
way, this has no effect on the image quality (see Eq. ([5])). 
In Figure[3]we compare GI of the Weiss domains with tra¬ 
ditional Faraday imaging (TFI); for the chosen angle of 
the polarizer P2 we observe the intensity distribution 
on the CCD area illuminated by a-beam, before summing 
pixels values for the bucket detector. The images of the 
magnetic domains of our sample obtained with GI tech¬ 
nique (a) and with TFI (b), averaging over K = 46904 
frames, include a (400 x 560)/rTO^ area characterized by 
equally distributed transmittances, such that i?+ = R- 
= R. Each pixel of the ghost image corresponds to the 
value of the correlation Cx^ between Na and the Tth 
pixel measurement Weiss domains shape and po¬ 
sition are evident in GI reconstruction, showing a good 
agreement with TFI. Differences between the two images 
are mostly related to a reduction of the resolution in GI. 
In TFI, the resolution of the image corresponds to the 
physical pixel dimension (8 /im) while for GI we evalu¬ 
ated FWHMc(^-) ~ 3 pixels. Figured shows the SNR of 
the ghost image as a function of the number of speckles 
R collected in areas with transmittance T'a_+, showing a 
clear agreement with the theoretical prediction. Red line 
refers to a theoretical model of SNR similar to Eq. |8]but 
considering also a background contribution. In fact, in 
our case, since we estimate Ta^+ = 0.028, Ta^- = 0.011 
and our pseudo-thermal light with M = 1 provides a 
number of photons of A = 1210, the detected mean inten¬ 
sity are quite low: = 34 for G Sa,+ and rix„ = 13 

for Xq G Sa^-. Therefore, background variance contribu¬ 
tion {Vback = 151) is not negligible in GI reconstruction. 

Conclusion. In this paper we proposed the application 



R 

FIG. 4. SNR curve as a function of the number of resolution 
cell collected in the area with transmittance Ta,+ . Red line 
correspond to the theoretical prediction. 


of the GI technique to magneto-optical Faraday (or Kerr) 
microscopy. GI is a well established imaging technique 
in the realm of quantum optics, that allows retrieving 
the image of an object from an optical beam it has never 
interacted with. This is achieved by exploiting the cor¬ 
relation with its conjugated beam interacting with the 
object and observed with a bucket detector. Due to its 
flexibility, this technique appears particularly promising 
if magnetic domains imaging has to be performed in hos¬ 
tile environment (very high magnetic field, sub-K regime, 
limited optical access to the sample), where the possi¬ 
bility of collecting the light after the interaction with 
the sample in a minimal amount of space, even without 
spatial resolution, e. g. by means of a single optical 
fiber, allows to overcome the limitations of the experi¬ 
mental setup. To prove the validity of this technique in 
the field of MOI we performed a first proof-of-principle 
experiment exploiting the Faraday effect for imaging the 
domains of a YIG sample. The techniques is extremely 
flexible and both Kerr or Faraday configurations can be 
exploited. We achieved a resolution of 24 /rm by means 
of a purpose developed microscope but, with a proper de¬ 
sign, the resolution of the ghost image can be increased at 
the level of traditional imaging. Our results demonstrate 
that this technique is able to provide high quality images 
of magnetic domains. In future, a system devoted to GI 
of magnetic domains to be used in hostile environment 
will be developed. This will be based on the configura¬ 
tion known as ”computational-ghost-imaging” [s^, 1^ . 
using a fibre-bundle to bring the multi-speckle-like light 
on the object, and a pigtailed photodiode with a multi- 
mode fiber as the bucket detector. 
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